First-principles calculations of phase transition, elasticity, and thermodynamic 

properties for TiZr alloy 
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Structural transformation, pressure dependent elasticity behaviors, phonon, and thermodynamic 
properties of the equiatomic TiZr alloy are investigated by using first-principles density-functional 
theory. Our calculated lattice parameters and equation of state for a and uj phases as well as 
the phase transition sequence of a— >uu— >/3 are consistent well with experiments. Elastic constants 
of a and w phases indicate that they are mechanically stable. For cubic f3 phase, however, it 
is mechanically unstable at zero pressure and the critical pressure for its mechanical stability is 
predicted to equal to 2.19 GPa. We find that the moduli, elastic sound velocities, and Debye 
temperature all increase with pressure for three phases of TiZr alloy. The relatively large B/G 
values illustrate that the TiZr alloy is rather ductile and its ductility is more predominant than 
that of element Zr, especially in /3 phase. Elastic wave velocities and Debye temperature have 
abrupt increase behaviors upon the a— >ui transition at around 10 GPa and exhibit abrupt decrease 
feature upon the uj— Yf3 transition at higher pressure. Through Mulliken population analysis, we 
illustrate that the increase of the ci-band occupancy will stabilize the cubic ft phase. Phonon 
dispersions for three phases of TiZr alloy are firstly presented and the j3 phase phonons clearly 
indicate its dynamically unstable nature under ambient condition. Thermodynamics of Gibbs free 
energy, entropy, and heat capacity are obtained by quasiharmonic approximation and Debye model. 

PACS numbers: 62.20.de, 64.70. kd, 61.50.Ks, 74.70.Ad 
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I. INTRODUCTION 



Group IV transition metals such as titanium, zirco- 
nium, and hafnium as well as their alloys have been 
investigated extensively by experiments and theoreti- 
cal calculations since their particular applications in the 
aerospace, nuclear, and chemical industries [lHjj. Their 
properties of low thermal neutron absorption cross sec- 
tion, adequate strength and ductility, good corrosion re- 
sistance, long-term dimensional stability in an irradia- 
tion environment, excellent compatibility with the fuel 
and coolant, good thermal conductivity and adequate 
resistance to fracture project them to the front of ma- 
terial applications for both fuel element cans and in-core 
structural components in water-cooled nuclear reactors 
[l| . Specifically, TiZr alloy is traditionally used for high- 
pressure neutron diffraction and has similar machining 
properties to the ZrNb alloy used in nuclear-reactor tub- 
ing. Therefore, investigations of their structural and elas- 
tic properties, which relate to various fundamental solid- 
state properties such as interatomic potentials, equation 
of state, phonon spectra, and thermodynamic properties, 
have prominent meaning in both scientific and technolog- 
ical fields. The scientific interest in Group IV transition 
metals and their alloys origins from the fact that they 
have a narrow d-band in the midst of a broad sp-band. 
The rf-band occupancy is crucial for the structural stabil- 
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ity of these materials and a pressure-induced s-d electron 
transfer is the driving force behind their structural and 
electronic transitions [__, [_} . 

At ambient condition, both Ti and Zr crystallize in 
hexagonal closed-packed (hep) structure (a phase). At 
high temperature of 1155 and 1135 K, they transform 
into the body-centered cubic (bec) f3 phase, respec- 
tively. At room temperature and under compression, a- 
Zr transforms to another hexagonal structure uj phase 
at about 2—7 GPa f8rlll|. Under further high pressure 
of 30—35 GPa, the uj— >/3 phase transition has been ob- 
served [13, [l_H___l ■ The whole a— >w— transition series 
for Zr have been reproduced in theoretical investigations 
[l6l - [l9j . However, the experimental established room 
temperature transition sequence of Ti is a— >uj— >"f— >5 
[lEl [H UJ. The /3 phase of Ti metal has not been ob- 
served up to 216 GPa f2fj|| . R ecent first-principles study 
performed by Mei et al. [22| found that the S-Ti is not 
stable under hydrostatic compression and the K phase 
transition sequence of Ti is a— >us— >"f— >/3. The absence of 
the high-pressure /3 phase for Ti in experiments was at- 
tributed to the possible nonhydrostatic stress which dis- 
torts the /3 phase [HI]. 

For Ti— Zr alloys, transition sequence of a— >uj^>-/3 upon 
compression has been reported by Bashkin et al. [24 - 
|26| |. Using the differential thermal analysis (DTA) and 
calorimetric technique, they extensively studied the ef- 
fects of temperature and pressure on the a^-ui and uiT^-fj 
transformations in the equiatomic TiZr alloy at temper- 
atures up to 1023 K, and pressures up to 7 GPa. [24[ 
At atmospheric pressure, the transition tempera- 
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ture was measured to be 852 K upon cooling. Under 
pressure, the transition temperature decreases down to 
the triple a— (3— uj equilibrium point (P tr =4.9±0.3 GPa, 
T tr =733±30 K). At pressures above the triple point, the 
slope of the ui~ ft equilibrium boundary is positive at 
pressures up to 7 GPa. At room temperature, the ol^-uj 
transition occurs at around 5 GPa. Cooling of the f3 
phase in the pressure range 2.8—4.8 GPa can form a two- 
phase mixture of a stable a and a metastable uj phase. 
Another group [27| reported that the equilibrium a—ui 
boundary is situated on the P — T diagram at 6.6 GPa. 
Afterward, Bashkin et al. [Ill [26[ employed the energy- 
dispersive X-ray diffraction (EDXD) with synchrotron 
radiation to investigate the structural transitions of the 
TiZr alloy. They found an increase in the ui— >/3 tran- 
sition pressure from about 30 to 43—57 GPa when the 
titanium content in the alloys increases from to 50 
at. %. The pressure dependant behaviors of supercon- 
ductivity was also presented. For equiatomic TiZr al- 
loy, they found that the a phase remains the sole stable 
phase under quasi-hydrostatic pressure up to 12.2 GPa 
26]. Only from 15.5 GPa on, the u> phase becomes dom- 
inant. Recent in situ high-temperature high-pressure 
angle-dispersive synchrotron radiation diffraction exper- 
iment revealed that under pressure the equiatomic TiZr 
alloy occurs an uj— ¥uj' isostructural transition at high- 
pressure domain (28| . 

However, to date only one theoretical study has fo- 
cused on TiZr alloy [231, where they calculated the 
electron structure and total energy by the scalar rel- 
ativistic full-potential linearized augmented-planewave 
(FPLAPW) method with a generalized gradient approx- 
imation (GGA) potential. The P — T phase diagram and 
tendency toward decomposition in equiatomic TiZr al- 
loy were calculated within the electron density functional 
theory (DFT) and the Debye-Griineisen model. Their 
calculations showed that the uj phase is stable at atmo- 
spheric pressure and moderate (up to 610 K) tempera- 
tures, and hence no a— ¥ui transition occurs at room tem- 
perature. They attributed the discrepancy between ex- 
periment and theory to the experimental samples, within 
which defects may exist. In present work, we have care- 
fully calculated the structural, phase transition, elastical, 
phonon vibrational, and thermodynamic properties for 
equiatomic TiZr alloy by employing the first-principles 
total energy calculations. Elastic constants, elastic prop- 
erties, and phonon dispersions of three experimentally 
observed phases are firstly presented. Our main motiva- 
tion is to give out detail pressure behaviors of the elastic- 
ity as well as the thermodynamic properties to support 
the practical application of TiZr alloy in nuclear technol- 
ogy 



II. COMPUTATIONAL METHODS 

First-principles DFT calculations on the basis of the 
frozen-core projected augmented wave (PAW) method 



of Blochl [30j are performed within the Vienna ab ini- 
tio simulation package (VASP) [3111 . where the Perdew, 
Burke, and Ernzerhof (PBE) |32| form of the GGA is 
employed to describe electron exchange and correlation. 
To obtain accurate total energy and stress tensor, a cut- 
off energy of 500 eV is used for the plane-wave set. The 
T-centered k point-meshes in the full wedge of the Bril- 
louin zone (BZ) are sampled by 18x18x16, 16x16x9, 
and 18x18x18 grids according to the Monkhorst-Pack 
(MP) 33} for a (two atoms cell), uj (six atoms 1x1x2 
supercell), and j3 (two atoms cell) phases TiZr alloy, 
respectively. Full geometry optimization at each vol- 
ume is considered to be completed when all atoms are 
fully relaxed until the Hellmann-Feynman forces becom- 
ing less than 0.001 eV/A. The Ti 3s 2 3p 6 3d 3 4s 1 and the 
Zr 4s 2 4p 6 4<i 3 5s 1 orbitals are explicitly included as va- 
lence electrons. 



III. RESULTS 

A. Structure 

In calculations of a phase TiZr alloy, we use one lamella 
structure, which has the same unit cell as the conven- 
tional hep unit cell but the Ti and Zr atoms are located 
on the layer along the c axis alternately. For j3 phase, 
Ti atoms are fixed at corner and Zr atoms are located at 
center of the bec unit cells. For uj TiZr, two types (A- 
type and B-type) of 1x1x2 supercells are constructed 
(Fig. 1). The theoretical equilibrium structural parame- 
ters, bulk modulus B, and pressure derivative of the bulk 
modulus S obtained by fitting the energy-volume data 
in the third-order Birch-Murnaghan equation of state 
(EOS) [HI are tabulated in Table 1. Experimental values 
from Ref. [26| and other theoretical results from Refs. [l[ 
and [2t| are also presented for comparison. Note that we 
only list in Table 1 the results of B-type structure for uj 
phase. Our calculated equilibrium structural parameters 
for A-type uj phase are ao=4.829 A and co=3.013 A. Its 
bulk modulus B and pressure derivative of the bulk mod- 
ulus B 1 are equal to 100.6 GPa and 3.63, respectively. All 
these values are almost equal to that of B-type structure. 
However, the total energy of B-type structure is smaller 
by 0.0769 eV per formula unit than that of A-type struc- 
ture. Using A-type structure for uj phase will lead to 
incorrect results in calculations of phase transition and 
mechanical properties. Thus, the results of A-type ui 
phase are not presented in the following and unless oth- 
erwise stated results of ui phase refer to that of B-type 
structure. 

From Table 1, one can find excellent coincidence be- 
tween our calculated values and the corresponding ex- 
perimental results for a and uj structural parameters. 
But our calculated bulk moduli B for these two phases 
are largely smaller than that from experiment. Although 
we can not simply attribute these results to the under- 
binding effect of the GGA approach, we find that our 
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TABLE I: Calculated lattice parameters (a or c), bulk modulus (B), pressure derivative of the bulk modulus (B ), and elastic 
constants of a, cj, and f3 TiZr at different pressures. For comparison, experimental values and other theoretical results are also 
listed. 
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FIG. 2: (Color online) Calculated enthalpy differences of a 
and uj phases with respect to /3 phase as a function of pressure. 



calculated bulk modulus B for a TiZr alloy lies within 
the range of experimental values between Ti 102 GPa |i| 
and Zr 92 GPa II]. As for (3 phase, our calculated equi- 
librium crystal constant is in good agreement with recent 
theoretical calculation (29|, where a value of 3.417 A was 
given out. The discrepancy between our results and the 
experimental value (3.417 A at 57 GPa) (2(| mainly due 
to the temperature effect. 



B. Phase transition at K 

At K, the Gibbs free energy is equal to the enthalpy 
H. After calculation, we can plot in Fig. 2 the relative 
enthalpies of the a and the U) phases with respect to 
the B phase as a function of pressure. The crossing be- 
tween the a and oj enthalpy curves readily gives phase 
transition pressure of —13.8 GPa, which indicates that 
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TABLE II: Calculated transition pressure for TiZr. For com- 
parison, other theoretical works and available experimental 
data for TiZr, Ti, and Zr are listed. 
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FIG. 3: (Color online) Calculated compression curves of TiZr 
alloy compared with experimental measurements from Ref. 
[25]'. 



at ambient pressure the u> phase is more stable than the 
a phase. This fact is in disagreement with experiment 
[24 126H27I , but coincides with recent theoretical calcula- 
tions [29t] Trubitsin et al. [29| attributed the discrepancy 
between experiment and theory to the experimental im- 
pure samples. They presented detailed explanation for 
this discrepancy. In our opinion, the imperfect crystals 
used in experiments may responsible mainly for this dis- 
crepancy. After all, the metastable u> phase was obtained 
at atmospheric pressure through cooling of the j3 phase 
under a pressure of 6 GPa with subsequent unloading at 
room temperature. One other reason may from the tem- 
perature contribution. In Table 2, we list the transition 
pressures for TiZr alloy as well as pure Ti and Zr metals 
from experiments and theoretical calculations. For the 
stable phase of metals Ti and Zr at ambient pressure, 
there exist debates in theoretical studies [ljl [lTj, 
Irij although experiments 0, [HO, H(| have reported the 
most stable phase to be a phase. For metal Zr, two 
DFT works [3a, |36j using full-potential linear muffin-tin 
orbital (FPLMTO)-GGA have illustrated the a^oj tran- 
sition for Zr at 4.8 and 3.3 GPa, respectively. However, 
two other DFT works usin § FPLMTO-GGA and 

FPLAPW-LDA methods indicated that the most sta- 
ble phase at K under ambient pressure is uj, not a 
phase. Two recent DFT-PBE works [H 03 also sup- 
port this conclusion. For metal Ti, the debates also ex- 
ist in DFT calculations [22|, |39| as indicated in Table II. 
By adding the phonon contribution to the Gibbs energy, 
recent DFT-PBE 0, E3 studies explicitly showed that 
the transition of a— >ui for Ti (Zr) occurs at 1.8 (1.7) GPa 
at room temperature. Thus, the disagreement between 
some theoretical works and experiments for elemental Ti 
and Zr mainly originates from temperature effect. The 
entropy from the thermal population of phonon states 
will stabilize the a phase at room temperature (22| . In 
our present work, we should also consider the tempera- 
ture effect on the a^-uj transition pressure for TiZr alloy 
since the same scheme used with respect to our previ- 




Pressue (GPa) 

FIG. 4: (Color online) Calculated lattice parameters in com- 
parison with the experimental data of Bashkin et al. [2(| for 
a and to phases of TiZr alloy. Our calculated results of uj' 
phase are also presented. 



ous study of Zr (T^- In addition, Trubitsin et al. [29| 
observed that a TiZr alloy is stable in the temperature 
range 600 K<T<900 K after including lattice energy and 
entropy in the Debye model. Therefore, the disagreement 
of a— >uj transition pressure between our calculation and 
experiment for TiZr alloy can also be attributed mainly 
to the experimental impure samples and partially to the 
temperature effect. As for the uj— >j3 transition, our calcu- 
lated transition pressure (33.9 GPa) is in good agreement 
with experiment [261 ]. 

To further analyze the pressure behavior of TiZr al- 
loy, we plot in Fig. 3 the compression curves of a, u>, 
and (3 phases. For comparison, the experimental data 
from Ref. (HI are also shown in the figure. Clearly, our 
calculated P— V equation of states are well consistent 
with the experimental measurement. The slight small- 
ness of our calculated volumes compared to that of the 
experimental data can be attributed to the temperature 
effect. The experiment was performed at room temper- 
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FIG. 5: Calculated moduli for TiZr alloy as a function of 
pressure. Solid lines are guides to the eyes. 



ature, while our data are valid only at K. In Fig. 4, 
we compare the calculated lattice parameters of a and 
uj phases with experimental results [26] . Table 1 also 
lists the optimized structural parameters of TiZr alloy at 
some different pressures. In the whole pressure domain 
of 0—60 GPa, the predicted lattice parameters for both 
a and uj phases compare well with the experimental val- 
ues, which supplies the safeguard for our following study 
of mechanical and elastic properties of TiZr alloy under 
pressure. 



C. Elasticity 

Elastic constants can measure the resistance and me- 
chanical features of crystal to external stress or pressure, 
thus describing the stability of crystals against elastic 
deformation. In present work, elastic constants for TiZr 
alloy in a, uj, and /? phases at different pressures are 
calculated through applying stress tensors with various 
small strains onto the optimized structures. The strain 
amplitude S is varied in steps of 0.006 from S=— 0.036 
to 0.036. Results are presented in Table 1. Obviously, 
the mechanical stability of a and uj TiZr at GPa and 
at some typical finite pressures can be predicted from 
the elastic constants data. But the elastic constants of /3 
phase at GPa illustrate that the cubic phase is mechan- 
ically unstable. We notice that the /3 phase metals Ti 
and Zr are also mechanically unstable at GPa [TH, |43[ . 
Along with the increase of pressure from GPa to 60 
GPa, the value of Cn — C12 increases near linearly from 
-5.3 GPa to 146.3 GPa for TiZr. Fitting the curve (not 
shown) of the pressure behavior of Cn — C12 by poly- 
nomial function, we find that the value of Cn — C12 be- 
comes positive at P >2.19 GPa. In fact, at P=2.5 GPa 
our calculated Cn,Ci2, and C44 equal to 102.9, 101.2, 
and 25.6 GPa, respectively, which explicitly indicate the 
elastically stable of bcc TiZr under this pressure. 

After obtaining elastic constants, the polycrystalline 




10 20 30 40 50 60 
Pressure (GPa) 



FIG. 6: Transverse and longitudinal sound velocities (a) and 
Debye temperature (b) for TiZr alloy as a function of pressure. 
Solid lines are guides to the eyes. 



bulk modulus B and shear modulus G are calculated 
from the Voigt-Reuss-Hill (VRH) approximations [iij . 
Results of TiZr alloy in a, uj, and /3 phases are shown 
in Fig. 5. The deduced bulk moduli from VRH approx- 
imations for all three phases at GPa turn out to be 
very close to that obtained from the EOS fitting, which 
indicates that our calculations are consistent and reli- 
able. Along with the increasing of pressure, both bulk 
modulus B and shear modulus G increase near linearly 
for all three phases. The increasing rates of B for all 
three phases are apparently larger than that of G. It is 
well known that the shear modulus G represents the re- 
sistance to plastic deformation, while the bulk modulus 
B can represent the resistance to fracture. A high (low) 
B /G value is responsible for the ductility (brittleness) of 
polycrystalline materials. The critical value to separate 
ductile and brittle materials is about 1.75 [45] . Using the 
calculated values of bulk modulus B and shear modulus 
G for TiZr alloy, the B/G values for a phase increase 
from 2.80 to 3.04 under pressure from GPa to 10 GPa, 
for uj phase increase from 2.08 to 2.57 upon compression 
from GPa to 40 GPa, and for f3 phase decrease from 
11.57 to 3.01 under pressure from 2.5 GPa to 60 GPa. 
These results indicate that transit to uj phase will lead to 
less ductility and further transit to (3 phase will result in 
more predominance of ductility. Although the large B/G 
values for f3 TiZr at pressure range of 0—30 GPa can not 
be achieved in experiment at room temperature, the B/G 
value of about 3.5 for (3 phase at 40 GPa indicates that 
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FIG. 7: Number of electrons on s, p, and d orbitals for TiZr 
alloy as a function of pressure, relative to their values at 
pressure. The numbers of s electrons (N s ), p electrons (N p ), 
and d electrons (Nd) at P=0 GPa for Ti atom are 2.469, 
6.417, and 3.012, respectively, for Zr atom are 2.565, 6.413, 
and 3.012, respectively. 

the ft phase TiZr alloy possess the biggest ductility. For 
a Zr, the B/G value of 2.63 at ambient condition can 
be derived from the elastic data in Ref. [46j. Therefore, 
TiZr alloy is rather ductile and its ductility is more pre- 
dominant than that of element Zr. We hope that our 
calculated elastic constants and elastic moduli can be il- 
lustrative in the realistic application of the mechanical 
data for TiZr alloy. 

Elastic wave velocities provide important information 
about the behavior of materials before phase transi- 
tion. The transverse and longitudinal elastic wave ve- 
locities of the polycrystalline materials can be calculated 
through relations of v t = y/G/p (p is the density) and 
v\ = y/ (3B + AG)/3p, respectively. The average wave 
velocity in the polycrystalline materials is approximately 
given as v m = [(l/3)(2/v| + 1/if )] -1 / 3 . Using the re- 
lation 9]j = (/i/fcs)(3n/47rr2) 1 / 3 z; m , the Debye temper- 
ature (0d) can be obtained. The calculated results of 
transverse and longitudinal elastic wave velocity are plot- 
ted in Fig. 6(a) and the Debye temperature in Fig. 6(b). 
Increasing behaviors of elastic wave velocities and 9d 
under pressure are obvious. Transiting from a to uj at 
around 10 GPa (experimental transition pressure), TiZr 
alloy exhibit abrupt increasing feature for both transverse 
and longitudinal elastic wave velocities: from 5.43 to 5.76 
km/s for v t and from 2.60 to 3.07 km/s for v\. In study of 
element Zr, Liu et al. [46j observed similar behaviors of 
elastic wave velocities by using ultrasonic interferometry 
in conjunction with synchrotron x-ray radiation. From 
Fig. 6(a), we find that the elastic wave velocities will 
decrease upon the uj to (3 transition at high pressure. 
The Debye temperature also shows this kind of increase 
first (at a to uj transition) and decrease later (at oj to /3 
transition) behavior upon compression. 

For discussion of pressure induced s-d electron trans- 
fer for TiZr alloy, we have performed Mulliken population 



analysis [47j of its /? phase and then plot in Fig. 7 the 
variation of the number of electrons on s, p, and d or- 
bitals with increasing pressure. Clearly, only the Nd of 
Ti atom increases with pressure. Arithmetic average of 
Nd for Ti and Zr atoms also possesses obvious pressure- 
dependent increasing behavior. Increase of d-band occu- 
pancy will stabilize j3 phase of TiZr alloy under pressure. 
However, we find that the sp electrons for Ti atom and 
the spd electrons for Zr atom all decrease upon compres- 
sion up to 143 GPa. No increasing behavior of N s has 
been observed. This fact is different from Ti [43| and 
Zr [48| , where it was found that the N s in both of these 
two elemental metals start to increase when the pressure 
exceeds about 70 and 100 GPa, respectively. Thus, the 
s-d electron transfer behaviors for TiZr alloy are different 
from its archetype metals. 



D. Phonon dispersion 

To our knowledge, although plenty of attentions have 
been paid on the phonon dispersions, thermodynamics, 
and P - T phase diagrams of Ti 0, EH and Zr pJ-Ql 
[50( | . no experimental and theoretical phonon dispersion 
results have been published for equiatomic TiZr alloy. 
In this subsection and the following subsection, we want 
to present our calculated results of phonon dispersion 
curves, thermodynamics, and P — T phase diagram for 
TiZr alloy. 

We use the supercell approach [5lJ and the small dis- 
placement method as implemented in the FROPHO code 
[52j to calculate the phonon curves in the BZ and the 
corresponding phonon density of states (DOS). In the 
interpolation of the force constants for the phonon dis- 
persion curve calculations, 5x5x5, 3x3x1, and 5x5x5 
fc-point meshes are used for a 3x3x3, uj 3x3x3, and (3 
4x4x4 supercells, respectively. The forces induced by 
small displacements are calculated by VASP. 

The calculated phonon dispersion curves and the 
phonon DOS of a, uj, and (3 TiZr alloy at GPa are 
displayed in Fig. 8. Overall speaking, behaviors of the 
phonons in the BZ for a and uj TiZr alloy are similar 
to that of metal Ti [H| and Zr [l6j]. Same with Zr, 
the uj phonons for TiZr alloy are stiffer along the c axis 
than in the basal plane due to the low c/a ratio. The 
a phonons along [001] direction for TiZr alloy are softer 
than that of metal Ti and Zr. Compared with (3 Ti and 
f3 Zr, phonons in the BZ for (3 phase of TiZr alloy ex- 
hibit different behaviors. The unstable phonon branch 
along [110] direction in the two metals transfer to show 
stable and stiff feature in the alloy. On the contrary, the 
phonons along [001] direction are stable for metals but 
unstable for the alloy. These informations are important 
in analyzing phase transformation. The stable nature 
of phonon branch along [110] direction for (3 TiZr alloy 
may responsible for the difficulties of a— ¥f3 transition un- 
der high temperature and the uj— >(3 transition under high 
pressure as indicated in experiments [H], HH, [27j ■ 
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FIG. 8: Phonon dispersion curves and phonon DOS of a, uu, 
and P TiZr alloy. 



E. Thermodynamic properties and P — T phase 
diagram 

Thermodynamic properties of equiatomic TiZr alloy 
can be determined by phonon calculation using the quasi- 
harmonic approximation (QHA) [53l. l54j or by the quasi- 
harmonic Debye model [55]. Within these two models, 
the Gibbs free energy G(T,P) is written as 



G(T, P) = F(T, V) + PV. 



(1) 



Here, F(T,V) is the Helmholtz free energy at tempera- 
ture T and volume V and can be expressed as 

F(T,V) = E(V) + F vib (T,V) + F el {T,V), (2) 

where E(V) is the ground-state total energy, F v n,(T,V) 
is the vibrational energy of the lattice ions and F e i (T, V) 
is the thermal electronic contribution. 

Under QHA, the F v n,(T, V) can be calculated by 



F vib (T,V) = k B T 







/ sM In 


2 sinh — 


/o 





dui, 



where uj represents the phonon frequencies and g(uj) is 
the phonon DOS. This formula requests positive results 
of phonon DOS. So it is not suitable for dynamically 
unstable phases. Instead, the vibration energy for phases 
with imaginary phonon frequencies can be estimated by 
the Debye model 



F V ib(T, V) = h B e D +k B T 



3 In 



(3) 



(4) 

where ^k B 8o is zero-point energy due to lattice ion vi- 
bration at K and D(6 B /T) the Debye integral written 

as D(0 D /T) = 3/(9 D /T) 3 J® d/T x 3 /{e x - l)dx. Detailed 
computation scheme of Debye model please see Ref. [55[ . 

F e i in Eq. (2) can be obtained from the energy and 
entropy contributions, i.e., E e i e — TS e i e . The electronic 
entropy S e i e is of the form 

S e ie(V,T) = -k B J n(e, V) [/ In /+(l-/)ln(l -/)]&, 

(5) 

where n(e) is electronic DOS and / is the Fermi-Dirac 
distribution. The energy E e i e due to the electron excita- 
tions takes the following form: 

E ele (V, T) = Jn(e, V)fede - J " n(e, V)ede, (6) 

where Ef is the Fermi energy. 

Using QHA and Debye model, we calculated the Gibbs 
free energy (G), entropy (S*), and specific heat at con- 
stant volume (CV) for a, uj, and /? phases of TiZr al- 
loy at their equilibrium volumes. As shown in Fig. 9, 
the Gibbs free energy and entropy calculated by QHA 
are almost identical to those obtained by employing the 
Debye model for a phase, while the differences between 
these two schemes for uj phase is slightly large, especially 
for high-temperature domain. Although the phonon cal- 
culation is more reliable in general analysis, some non- 
harmonic terms, which is important for u> phase, possi- 
bly having been ignored in the QHA method. From Fig. 
9(a), cross between a and uj curves by Debye calculations 
readily gives phase transition temperature of 680 K. A 
phase transition temperature of 940 K for a— >(3 also can 
be seen. So, we can use Debye model to obtain the P — T 
phase diagram by changing the pressure P (see Fig. 10). 
However, there is no cross between a and uj curves us- 
ing phonon calculations. This need more works to clarify. 
Figure 9(b) shows the entropy results. While the entropy 
of lu phase is lower than a phase, the results of /? phase 
are higher than a phase in all considered temperature 
domain. The specific heat at constant volume (Gy) can 
be directly calculated by T(|^)y. Results are presented 
in Fig. 9(c), using which one can estimate the Debye 
temperature. Actually, our phonon calculations indicate 
that the zero-point energies for a and uj phases at P=0 
GPa are 27.02 and 28.00 meV, respectively. So, the De- 
bye temperature values from phonon calculations equal 
to 278.6 and 288.9 K, respectively. The Debye model 
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FIG. 9: Temperature dependences of (a) Gibbs free energy, (b) entropy, and (c) specific heat at constant volume (Cy) for a, 
u, and j3 phases of TiZr alloy. 
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FIG. 10: Calculated P -T phase diagram of TiZr alloy. For 
comparison, we also plot in figure the experimental data [24[ 
and the theoretical results from Ref. [29l ]. 

gives 6?£>=303.6 and 369.0 K for a and u phases, respec- 
tively. So, the Debye model gives closer values than the 
phonon calculations with respect to the Debye temper- 
ature values of 305.1 and 355.6 K computed by elastic 
constants (see Fig. 7). 

We present in Fig. 10 the P — T phase diagram of 
equiatomic TiZr alloy calculated by the Debye model. 
A triple point (P tr =2.5 GPa, T tr =760 K), comparable 
to experimental results (P tr =4.9±0.3 GPa, T tr =733±30 
K) [24( and previous theoretical values (Pt r =4.2 GPa, 
Ttr=720 K) [23], is predicted. The to— >a and a— >@ transi- 
tion temperatures of 680 and 940 K are close to previous 
Debye-Griineisen values (610 and 943 K, respectively) 
[2^ | . Detailed understanding of the P — T phase diagram 
please see Ref. [29l |. 

IV. CONCLUSION 

In summary, structure, phase transition, elasticity, 
phonon, and thermodynamics of the equiatomic TiZr al- 



loy have been studied by means of the first-principles 
DFT-PBE method. The calculated lattice parameters 
and equations of state for a and u) phases are in good 
agreement with experiments. We have found that the u> 
phase is most stable at GPa and the oj— >j3 transition 
was calculated to occur at around 33.9 GPa. The elas- 
tic constants of a phase TiZr alloy at P=Q GPa coincide 
well with previous calculations. Mechanical stability of a 
and uj phases were predicted and the (3 phase is mechani- 
cally unstable at zero pressure. Our elastic constants in- 
dicate that the /3 phase will become mechanically stable 
upon compression. Under pressure, both bulk modulus B 
and shear modulus G increase near linearly for all three 
phases. Increasing behaviors of elastic wave velocities 
and Debye temperature under pressure have also been 
predicted. It was found that elastic wave velocities and 
Debye temperature increase first (at a to w transition) 
and decrease later (at ui to j5 transition) abruptly upon 
compression. The Mulliken population analysis showed 
that either sp electrons of Ti atom or spd electrons of Zr 
atom transfer to the Ti-d orbital with increasing pres- 
sure. 

In phonon dispersion study, the dynamic stable nature 
for a and u> phases are observed. The unstable modes of 
j3 phase TiZr alloy along [001] and [111] directions have 
been shown and these imaginary vibration branches are 
different from elemental Ti and Zr. Based on phonon 
calculations and also the Debye model, thermodynamic 
properties of Gibbs free energy, entropy, and specific heat 
at constant volume as well as the P — T phase diagram 
for TiZr alloy have been presented. 
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